Lack of Self-affinity and Anomalous Roughening in Growth Processes 
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We contrast analytical results of a variety of growth models involving subdiffusion, thermal noise 
and quenched disorder with simulations of these models, concluding that the assumed self-affinity 
property is more an exception than a rule. In our two dimensional models, self-afRne surfaces may 
only appear when the roughness exponent is x = I /2 or x = I . A new scaling picture, which leads 
to more suitable ways of determining the scaling exponents, is proposed when lack of self-affinity 
exists. 



Interface roughening observed in many growth phenomena [1,2] can be explained as produced by the interplay 
between a surface tension, which tends to smooth the interfacial curvature, and some kind of noise coming either 
from thermal fluctuations or spatial inhomogeneities (i.e. quenched disorder). In the simplest model exhibiting 
kinetic roughening, the Edwards- Wilkinson (EW) model [3], both effects are involved by means of a diffusion term 
and an additive thermal noise. In order to introduce some symmetries existing in many real systems, a nonlinear 
term must be added to the EW model and this gives rise to the Kardar-Parisi-Zhang (KPZ) equation [4]. The two- 
dimensional KPZ model successfully describes a remarkable diversity of growth processes and is an excellent example 
of universality [1,2]. Simulations and analytical calculations indicate that the two-dimensional KPZ model leads to 
self-afEne interfaces with simple scaling laws. One of the most important quantities used to characterize the scaling 
of the interface is the global width: 



where L is the system size, < .. > denotes a spatial average over the whole system and {..} over realizations. For the 
early times regime, the interface width is expected to be of the form a(L,t) ~ t 13 where j3 is called the time exponent. 
At a characteristic saturation time t s (L) ~ L z , where z is the dynamic exponent, the horizontal correlation length 
reaches the system size and the width saturates, <r(L,t > t s ) ~ L x , with a roughness exponent x- This is the so-called 
Family- Vicsek scaling ansatz [1,2]. 

Thus, the roughness exponent is calculated from the scaling of the width in saturation. To measure x m an 
experiment or simulation, one has to consider systems with different sizes and wait for saturation at times t s (L) ~ L z . 
This way of determining x is n °t very convenient from a computational point of view due to the large system sizes in 
which simulations must be done. The situation in experiments is not better, because one should perform measurements 
in samples with different sizes. For that, in practice an alternative method is usually used. The width is evaluated 
over a window of size / <C L and one waits for saturation a shorten time t s (l) <C t s (L). This local width is believed 
to scale as the global width, cr(l,t) ~ t 13 for short times, and the roughness exponent is obtained from the scaling of 
<j(l, i) at times t l z : 



Equivalently, the height difference correlation function, G(l,t) = {((/}(/ + x,t) — h(x,t)) 2 )}, can be used to determine 
the scaling exponents because G(l,t) ~ a 2 (l,t). It is a general belief [1,2] that the roughness exponent xioc obtained 
in this way is equal to the exponent x that one gets from the global width. In fact, this is correct when the interface 
is self-afEne. 

In general, this picture has been extrapolated to other kinds of growth phenomena involving higher dimensions 
(KPZ in an arbitrary dimension) [1,2], subdiffusive behavior in epitaxial growth phenomena [5,2] and quenched noise 
in porous media [2,6-10]. In contrast to the two-dimensional KPZ case, now experiments and simulations show a great 
discrepancy. There are numerical evidences about an anomalous scaling behavior in growth processes in which the 
roughness exponent is larger than one (x > 1)- In numerical studies of models related to the molecular-beam epitaxy 
(MBE) growth, it has been found [11,12] that global roughness (measured by the interface width) gives different 
scaling exponents than local roughness (measured by either a height difference correlation function or the local width 
in a window of size /). Similar results have been numerically obtained in the case of the EW equation with quenched 
noise [8-10]. This equation has attracted much attention in recent years [2,6-10] because it is the simplest way of 
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describing driven interfaces in quenched disordered media. Surprisingly, it seems that local and global dynamic scaling 
behaviors are distinctly different for some growth processes. This anomalous scaling behavior is not predicted by the 
conventional dynamic scaling hypothesis of Family- Vicsek and it is not well understood. 

The most recent interpretation [8-11,13,14] of this anomalous scaling is that the conventional dynamic scaling fails 
in the case of roughness exponent larger than one , % > 1. In the case of MBE growth some other effects as the 
breaking of spectra [15] and the existence of a groove instability [11] have been also associated with anomalous scaling. 
However, very recently, Das Sarma et. al. [16] have studied MBE growth models with a global roughness x < 1 111 
which an anomalous scaling behavior appeared as well. So, it seems that a roughness x > 1 is not the only reason for 
anomalous scaling behavior. Moreover, the situation is further complicated due to fact that for most of these growth 
models there is not an analytical value of the roughness exponent. 

Our aim in this Letter is to show that anomalous scaling occurs not only whenever x > 1 but also in growth processes 
in which % < 1 different roughness exponents may appear at small and large length scales. As a consequence the 
interface is not self-affine and a new scaling picture emerges. Our conclusions are based on analytical solutions (and 
comparison with simulation) of growth models with subdiffusive behavior and/or quenched disorder. 

Models with uniform diffusion.- The inequality x 7^ Xioc has been numerically found in models without analytical 
solutions. In the following we are going to consider models with analytical solutions that will allow us to understand 
the existence of a local exponent Xioc- Let us start by studying models with diffusive or subdiffusive behavior and 
thermal noise governed by the equation: 

dt ~ [ ' dx 2 ™ +^ x ' t >' 

where n(x,i) is a Gaussian noise with correlation: 

(n(x,t)n(x' ,t')} = A a (x - x')6(t - t'), 

the correlator A a (x) is a normalized function decaying rapidly to zero over a finite distance a. The case m = 1 
corresponds to the EW model [3], for m = 2 we have the linear model [5,16,11,13] used in MBE growth and in general 
by increasing m more subdiffusive processes are obtained. 

Since these models are linear the scaling exponents can be easily calculated: x = (2m — l)/2, /? = (2m — l)/4m 
and z = 2m. We have numerically solved Eq.(3) for m = 1 and m = 2. In all cases there exists a good agreement 
between the analytically calculated x an d simulation results when the Eq.(l) was used. For m = 1 also the local 
method, Eq.(2), gives the correct exponent xioc = X = 1/2- However, for m = 2 there is a remarkable discrepancy 
between the local roughness exponent obtained in our simulation , 0.83±0.05, and the theoretical value x = 3/2. This 
discrepancy indicates that the interface is not self-affine and different roughness exponents exist at large and small 
length scales. The analytical treatment of (3) allows to understand the existence of this anomalous scaling. Consider 
the height difference correlation function G(l,t) under the l/L <C 1 condition: 

du ~ ^r^^-^ 11 * )A fl/tl /.(«) (4) 

-oo u 

A a (k) being the Fourier transform of the correlator A a (/). For early times, l/t 1 /* 1, we have G(l,t) ~ c(t)t 2 ^ , 
c(t) slowly varying and in the long time limit, \/t x l z < 1, G(l,t) ~ t 2 f> \\ /t 1 ' ' z )b(l / a) (where b(l/a) ~ constant) when 
m = 1 and G{l,t) ~ t 2 ^ (I /t 1 ! z ) 2 when m > 2. Then for small length scales, / <C t 1 ^ , we have the general scaling 

G(l,t) ~ t 2p {l/t 1 l z ) n ~ t 2l3 'l 2x '°*, (5) 

where /?* = [3 — n/2z, xioc = n/2, and the integer n = 1 or n = 2 depending on the model. Note that this scaling 
would occur whenever G(l,t) ~ t 2 ^ g(ljt 1 l z , a/t 1 ^), where the scaling function g(x,y) is analytical in x = and 
(/(O, 0) = 0. For true self-affine interfaces the exponent /?* = and the Family- Vicsek scaling is recovered. This only 
occurs for f3z = n/2, in such a way that the interface is self-affine whenever x = n /2 and in this case x = Xioc- Thus, 
Eq.(5) gives the self-affinity condition in the sense that whenever /?* ^ the interface is not self-affine. 

From Eq.(5) for m = 1 (case with n = 1), we obtain that G(l,t) becomes constant and independent of time, 
G{l,t ^> l z ) ~ /, as corresponds to the usual Family- Vicsek scaling picture. However, for m > 2 (n = 2 in this case) 
for small length scales, / <C t 1 ^ , we have G{l,t) ~ / 2 t 1-3 / 2m , and there is not saturation of G(l,t ^> l z ) at long times. 
Contrary to what is expected in the usual scaling, the height difference correlation (or equivalently the square local 
width) scales with time as t 213 * and a new temporal exponent appears, /?* = j3 — l/z = (2m — 3) /4m describing the 
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long time behavior of the local width. In particular, the linear model of MBE [5,16,11,13] (Eq. (3) for m = 2) exhibits 
an anomalous scaling in the sense that in the earliest time regime the local width (measured by the quantity G(l,t)) 
scales as t 2/3 ~ t 3 / 4 and does not saturate but G(l,t) ~ ft 213 * for / < t x l z < L (where /?* = j3 - l/z = 1/8). Our 
numerical results (see Fig. 1) for Eq.(3) with m = 2 are in excellent agreement with this analysis. The saturation 
regime in this case is dominated by the saturation of the largest scales / ~ L. At time t ~ L z saturation occurs in the 
whole system and G(l,t —> oo) ~ l 2 L 2/3 * z ~ l 2 L 2x ~ 2 for / <C L <C t 1 ! 7 ' . When saturation occurs, G(l,t —> oo) exhibits 
two different spatial scaling behaviors. For / <C L, G(l,t —> oo) scales with the window size as I 2 (which corresponds to 
a local roughness exponent xioc = 1) and for large scales, / ~ L, G(L,t — ► oo) ~ a 2 (L,t — ► oo) ~ L 2% . This crossover 
behavior leads to an effective exponent Xeff in between large and short length scales, that is the scaling behavior 
measured in simulations (in our simulation we found Xeff = 0.83 ±0.05). The local exponent Xeff apparently should 
take an intermediate value between the interpolated exponents xioc = 1 and x = 3/2 but it is never the case because 
G(l,t) is actually greater than the corresponding value if the interface were saturated and then the interpolation to 
largest scales occurs with decreasing slopes. As a consequence any determination of the local roughness exponent 
from Eq.(2) gives underestimated values of Xioc- 

Models with random diffusion.- As we have discussed above, models governed by Eq.(3) give anomalous scaling 
for m > 2. In this case the roughness is always x = (2m — l)/2 > 1 and xioc ~ 1 on short length scales. The 
question that we would like to discuss now is if the lack of self-affinity is due to a roughness exponent x > 1 or > on 
the contrary, anomalous scaling can exist even in growth models with x < 1- The most extended opinion [8-11,13,14] 
is that interfaces with x > 1 are not properly fractal objects and self-affinity is not a sure feature. But self-affinity 
for interfaces in which x < 1 is a simple hypothesis that, as we will show, in our models is only valid for x = 1/2 or 
1. This can be readily verified in a kind of growth model in which x takes a continuous range of values controlled by 
the intensity of disorder in the diffusion coefficient. This model can be used to describe a fluid interface advancing 
through stratified porous materials [17] with long range correlations in the growth direction. The growth model that 
we study here is given by 

D(x) > is the columnar diffusivity coefficient that controls the interconnectivity between neighbor channels. This 
equation can be viewed like (3) for m = 1 but with a quenched random diffusion term. A similar model was recently 
studied in Ref. [17]. 

It is known [17] that Eq.(6) can be considered as a model of diffusing particles with time-dependent sources and 
sinks in such a way that it is possible using the well known theory of transport through random media to take the 
usual stochastic description for D(x). So, D(x) is taken to be an uncorrelated random field and distributed according 
to a probability density P(D) = N a D~ a f c (D / D max ), where N a is a normalization constant and f c a cutoff function. 
The parameter a characterizes the intensity of disorder. Disorder is called weak for a = 0. Strong disorder occurs for 
< a < 1 because of the divergence of P(D) at origin. The model so formulated is easy to simulate numerically and 
allows analytical treatment. As it has been shown [17] , the inhomogeneous problem (6) is statistically equivalent to an 
homogeneous one with a time dependent diffusion coefficient given by the so-called Effective Medium Approximation 
condition [17]. With the above defined P(D) the equivalent diffusion length scales as lo(t) ~ t( 1-Q 0/( 2-Q and 
consequently the dynamical exponent z will be given by z = (2 — a)/(l — a). The exponents [3 and x ca n be easily 
calculated, obtaining [3 = 1/2(2 — a) and x = 1/2(1 — ot). Note that, as we have mentioned above, a continuous range 
of roughness exponents (from x = 1/2 to x = °°) are obtained by varying the intensity of disorder a. 

The height difference correlation function can be calculated [17] in the asymptotic time regime from (6) and is given 
by : 

!/•/•» „sv 

l/z 



G(l,t)~tVj_J U J o dv J d8 J o ^ (s _ s/ _ j(s _ s/)m2)(s/ _ j(s/)m2) (1-^ > a / t -(«), (7) 

where D(s) ~ s a /( 2-Q [ s the Laplace transform of the time-dependent diffusion coefficient. Since (6) has the same 
analytical behavior as (3) for m = 1 (with an equivalent time-dependent diffusion coefficient D(t)) it is easy to see 
that the scaling of G(l,t), Eq.(7), for small scales / <C t 1 ^ is also given by Eq.(5) with n = 1 (with the corresponding 
a-dependent exponents). 

For small length scales, / <C t 1 ^ <C L, from Eq.(5) we have that the square local width G(l,t) crosses over to a 
different temporal regime G(l,t ^> l z ) ~ t 2l3 *l where the time exponent is = j3 — l/2z = a/ 2(2 — a). Thus, the 
interface would be self-affine only in the case of weak disorder a = (i.e. the EW equation). For < a < 1, the 
interface is not a self-affine object and an anomalous scaling occurs. In Fig. 2 our numerical results for the local width 
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versus time are presented for three degrees of disorder. We can see that the exponent /?* predicted by our scaling is 
in very good agreement with simulation for all a (see Table I). 

In true saturation, at very long times, / <C L <C t 1 /* , the scaling behavior of G(l,t) can be obtained from (5) for 
n = 1 and we have G(l,t —> oo) ~ lL 2x ~ 1 . So, in the case of strong disorder a local roughness xioc = 1/2 is obtained. 
Certainly, as we can see in Table I, the crossover to the global roughness x = 1/P(2 — a)] > 1/2 for larger scales 
yields lower numerical values of the local roughness exponent. In all cases, included those with % < 1, the interface 
is not self-affine as we claimed above. From our models it is clear that the standard exponents x, ft, and z can not 
characterize completely the scaling behavior of a growing interface. For instance, note that the linear MBE model, 
Eq.(3) for m = 2, leads to the same critical exponents [3, a and z, as the model with random diffusion, Eq.(6), in the 
case of a disorder parameter a = 2/3, however both models have a very different local roughness exponent in such a 
way that they belong to distinct universality classes. 

In conclusion, we have dealt with a variety of growth models that allowed us to an analytical treatment. Analytical 
results have been compared to simulations focusing our interest on the scaling of the local width at small length 
scales, / <C L. We have concluded that, except for growth with x = n /2 (where n = 1, 2), interfaces in two dimensions 
are not self-affine and the Family- Vicsek dynamic scaling fails. In addition, we have seen that a new time exponent 
ft* = /3 — n/2z (where n = 1, 2 depending on the model) has to be introduced in order to get a complete description 
of the scaling behavior of the interface. Moreover, the roughness exponent Xeff , commonly obtained in a number of 
simulations (and experiments), is an effective exponent affected by a crossover from the local scaling exponent xioc to 
the global roughness exponent x- This effective roughness exponent does not give any information about the value 
of the scaling exponents x, ft and z, however, Xeff is important from a physical point of view to obtain the effective 
roughness measured in experiments and to distinguish the corresponding universality class. 

Although the new scaling picture, proposed in Eq.(5), has been obtained from the study of particular models, we 
believe that it can be applied to more general growth processes (in two dimensions). The scaling (5) is based on 
property of the scaling function, g, of being analytic at small length scales and this is expected to be verified for a 
broad class of growth problems. The existence of kinetic roughening satisfying (5) but with values of n ^ 1,2 must 
be investigated further since, in this case, xioc would be a new and interesting scaling exponent [18]. 
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TABLE I. Scaling exponents for the model with random diffusion, Eq.(6), for several degrees of disorder. Theoretical values 
of ft, a, and ft* are given in the text. The effective roughness exponent is always lower than the local roughness xioc = 1/2. 
Note that xioc remains roughly constant despite the great difference existing in \ 



a 


ft (exact) 


ft (simulation) 


X (exact) 


Xeff (simulation) 


ft* (theory) 


ft* (simulation) 


2/3 


3/8 


0.38 ± 0.05 


3/2 


0.35 ± 0.05 


1 / A 

1/4 


O.Z i ± 0.03 


1/2 


1/3 


0.34 ± 0.01 


1 


0.39 ± 0.04 


1/6 


0.18 ± 0.03 


1/3 


3/10 


0.31 ± 0.02 


3/4 


0.45 ± 0.05 


1/10 


0.12 ± 0.04 



FIG. 1. Square local width vs time for the linar MBE model, m = 2, calculated over a window size I = 5 (squares), I = 10 
(triangles), and I = 150 (crosses). The continuous line corresponds to the fit of the square global width, a 2 (L,t) ~ t 2/3 . The 
dashed line is the adjusted curve of the data for I = 5 and its slope 0.26 ± 0.02 corresponds to 2ft*. A system size of L = 1000 
was used and an average over 15 realizations of disorder was performed. 

FIG. 2. Square local width vs time for the random diffusion model, Eq.(6), averaged over a length scale I = 5 (squares), 
I = 10 (triangles), and I = 150 (crosses) for different degrees of disorder. A system size L = 1000 and 15 realizations of disorder 
were used. As in Fig.l, the continuous line slopes correspond to the values of 2ft obtained in simulation. The dashed line is 
the adjusted curve of the data for I = 5. A comparison with our scaling theory is given in Table I. 
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